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[57] ABSTRACT 

A method and system for the localization and charac- 
terization of the generators of human brain electromag- 
netic physiological activity includes a set of bioelectro- 
magnetic amplifiers, sensorial stimulators, and a com- 
puter based system for signal analog to digital conver- 
sion and recording. Sufficient statistics, including 
higher order statistical moments, for event related com- 
ponents are computed from the recorded signals, either 
in the time, frequency, or time-frequency domain, re- 
taining stationary, non-stationary, linear, and non-linear 
information. The localizations, orientations, activities, 
and connectivities of the generators are obtained by 
solving the inverse problem using sufficient statistics 
under anatomical, and functional constraints. Realistic 
head geometry and conductivity profiles are used to 
transform the measurements into infinite homogeneous 
medium measurements, through use of ananatomical 
deconvolution operator, thus simplifying optimally in- 
verse solution computations. Goodness of fit tests for 
the inverse solution are provided. Generator character- 
istics are visually displayed in the form of three and two 
dimensional head images, and optionally include proba- 
bility scaled images obtained by comparing estimated 
generator characteristics with those of a normal popula- 
tion sampled and stored in a normative data base. 

35 Claims, 3 Dramng Sheets 
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ing of the CNS. The signal recorded at each electrode is 
METHOD AND SYSTEM FOR THREE summarized by means of a set of descriptive parameters 

DIMENSIONAL TOMOGRAPHY OF ACTIVITY (DPs), based on stochastic process modeling. The DPs 
^J^^^DJ™ OF BRAIN **** mART reflect the normal and pathological functioning of the 
ELECTROMAGNETIC WAVES GENERATORS 5 CNS. Topographic maps based on the DPs are clini- 

** * /^/.t,^^^ _ , cally useful, and even more so when statistically com- 

BACKGROUND OF THE INVENTION pared t0 a normative data base. 

1. Field of the Invention However, this analysis method generates an exces- 
The present invention relates to electronic computer- sively large number of DPs, thus making quite difficult 

feed medical instruments and more particularly to the 10 the evaluation of a particular patient. Moreover, the 
localization and characterization of the generators of method does not attempt to localize the generators 
brainand heart electric and magnetic activity by a non- responsible for the measured DPs, thus limiting the 
invasive computerized method and system. clinical usefulness and contributing to the excessive 

2. Description of Related Art redundancy of the DPs due to volume conduction ef- 
The determination of the three dimensional localiza- 15 fects. Finally, EEG is limited to the study of second 

tion and of the temporal activity of the neuronal genera- order moments in the frequency domain, which means 
tors which give place to waveshapes, in an electroen- that the EEG has been implicitly assumed to be a 
cephalogram (EEG) and an magnetoencephalogram Gaussian stochastic process, despite evidence revealing 
(MEG) related to pathologies of the central nervous the non-linear nature of such signals, 
system (CNS), constitutes an important medical prob- 20 ln v s p at No 4(913160 a method for ^ reduction 
lem. Such knowledge can be helpful m producing more of the dimensionality of the DPs is proposed based on 

'? ^hT* * n ■ efSe ^T^l! 1 " 0 P™** 1 components (PCs) computation. This proce- 
K ^ treatments. A dure produces minimum ^ of ^ comb ination7of 

w,^k tf P ? t StUdy fOCU f jo FT ° m ori ^ DPs > ™ th optimum descriptive properties, 

The EEG and the MEG both have their common "L*ZT^r£ , ,T ^T*"* 0 * 
origin in the ionic currents produced by the cellular ™ Te > met ^ do ? P°\ take mto the non " 
elements (the neurons) composing the CNS. The total ^ear nature of the on^al sisals, 
current density vector field ^determined by the vecto- 30 ^Pfovement in the usefulness of QEEG has been 
rial additive combination of all of the elementary cur- «*«val b V means of biophysical models which take 
rents. The simultaneous activation of a large number of ^account the behavior of the electromagnetic fields 
such elements, together with an adequate geometrical P roduccd b J cm " rci J t sow** in a complex volume con- 
distribution, produces resulting electric potentials and ductor such M human head In this **** us - Pat- 
magnetic fields which can be measured outside the 35 N °s. 4,416,288; 4,753,246; and 4,736,751 propose proce- 
head. In the transformation process from total current dures for eliminating the distortion effects due to the 
density to measurable external fields, the effects of the volume conductor. However, they do not deal with the 
volume conductor properties of the different tissues spatial characterization of the generators, 
composing the head must be taken into account: brain, Several attempts have been made to fit equivalent 
meninges, cerebral spinal fluid, skull, and scalp. 40 dipoles *o measured fields in order to represent, albeit 

The resulting measured fields have the characteristics approximately, concentrated generators, either in the 

of a stochastic process, which can be described either in time domain or in the frequency domain. These proce- 

the frequency domain or in the temporal domain, as a dures are based oh the minimization of a certain distance 

function of the statistical moments. In the case of a criterion between the measurements and the theoretical 

Gaussian process, first and second order moments give 45 field values due to a current dipole inside a volume 

an exhaustive description. conductor model of the head. 

The neural elements which generate a given EEG or This type of procedure for source localization, based 

MEG component may be localized on a small cortical on ^ 1TSl order moment data, does not take into account 

area ("concentrated generator") or may, on the other the existence of diffuse generators, nor the existence of 

hand, be widely distributed in different parts of the 50 other sources of "spatial noise". Furthermore, a statisti- 

CNS ("diffuse generator"). The determination of the cal method for testing the goodness of fit of the source 

spatial distribution of the generators and of the multi- model is not provided. On the other hand, there is a 

variate statistical moments describing their interactions fundamental limit on the number of dipoles that can be 

is very important estimated, the maximum number being roughly equal to 

For a number of decades electric potential measure- 55 the number of electric or magnetic signals divided by 

ments of the CNS have been performed by means of six. 

electrodes placed on the scalp. Much experience has In French patent 2,622,990, several improvements are 
accumulated on the practical utility of the visual inspec- achieved by using frequency domain second order mo- 
tion of the EEG in the diagnostics and treatment of ment data, in the form of coherence matrices. An esti- 
patients with neuropsychiatry diseases. More recently, 60 mation method for the cross spectral spatial noise ma- 
brain magnetic fields have been measured (U.S. Pat. No. trix is proposed, under the assumption of interelectrode 
4,591,787), offering complementary information to that independence, the method thus being statistically equiv- 
obtained from the EEG. alent to the classical factor analysis model. The eigen- 
The current state of the art, as reflected in U.S. Pat. vectors of the common factor space are then used for 
Nos. 4,201,224; 4,408,616; and 4,417,592 is summarized 65 determining the concentrated generators (as many as 
as follows. Quantitative analysis of brain electric acti v- the number of common factors), 
hy by means of digital sjgnal processing methods However, empirical and theoretical evidence points 
(QEEG) allows an objective evaluation of the function- towards a diffuse generator model for spatial noise, 
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producing a structured cross spectral noise matrix for the corresponding neural tissue volume, for determin- 

EEG and MEG. This explains why the proposed noise ing the linear operator which transforms original EEG 

elimination method under the interelectrode indepen- and MEG measurements into equivalent infinite homo- 

dence assumption gives incorrect results. In such a case geneous medium measurements (anatomical deconvolu- 

computations based on coherence matrices are not justi- 5 tion); 

fied. Furthermore, dipole fitting methods applied to j) Using the parametric geometry for determining 
second order moment data or to eigenvector data are anatomical and functional constraints for the localiza- 
not equivalent Finally, interactions between genera- tions, orientations, activities, and connectivities of the 
tors, neither linear nor non-linear, are taken into ac- brain electromagnetic waves generators (generator con- 
count in the eigenvector dipole fitting approach. 10 straints); 

SUMMARY OF THE INVENTION k > Dfeital preprocessing of the EEG and MEG for 

. artifact and noise elimination, and for the separation of 

The objective of the present invention is a method EEG and MEG samples related to the fiducial markers, 

and system for the characterization of both concen- for obtaining event related components (ERCs); 

trated and diffuse generators of the EEG and MEG, t5 1) Statistical analysis of the ERCs for determining the 

based on all the statistical information available in these most adeqU ate numerical description of the spatio-tem- 

signals, in the form of statistical moments of all orders ^ properties m terms of suf r ici ent statistics; 

in the time or frequency domain. The mvention wdl m ) Computation of the activities and connectivities of 

allow the detection and carnation of the effect of the the ERCs generatorSj ^ on the static to the 

Affuse generators on the EEG and MEa Also, it will 20 Averse electromagnetic problem, under the above-men- 

allow the estimation of an increased number of concen- tQr ^ nst ^ ^ ^ tuffident statis . 

22g^ t0fiether Wlth their lmearand non - lm - tics for the ERCs transformed toln infinite homogene- 

In accordance with a first aspect of the invention OUS , medium * means °* ' thc anaton^cal deconvohition; 

there is provided a method for ihe three dimensional 25 n) in generator constraints do not altow a 

tomography of activity and connectivity of brain elec- ™« ue solutlo t n t0 * e ^ mmiber of 

tromagnetic waves generators, said method including: ™° generators should be decreased sufficiently to 

a) Attaching or approximating a set of electrodes and f low for the P™*** 'dentifiability of the mverse prob- 
magnetic 6ensors to the scalp of an experimental subject i « * L „ ^ „ , 
to detect brain electromagnetic physiological activity in 30 . °> evaluation of the goodness of fit of the 
the form of an electroencephalogram (EEG) and a averse solution, taking into account the existence of 
magnetoencephalogram (MEG), and measuring the ^lored spatial and temporal noise, and including statis- 
exact positions of the electrodes and sensors with re- tIcal hypotheses testing on the absence of activity and 
spec! to a reference coordinate system determined by connectivity of the ERCs generators; 

certain anatomical landmarks of the subject's head; 35 P) Optional computation of multivariate distances 

b) Amplifying the said electromagnetic signals de- between ERCs generators characteristics (localizations, 
tected at each electrode and sensor; orientations, activities, and connectivities) of the experi- 

c) Obtaining on-line digital spatio-temporal signals, mental subject and of a normal population as deter- 
consisting of the EEG and MEG, by connecting ana- mined from a normative data-base; 

log-digital converters to each amplifier, and digitizing 40 <0 Visual of th« dimensional and two di- 

al! data as it is gathered, under the control of a central mensional images corresponding to the localizations, 
experimental program; orientations, activities, and connectivities of the ERCs 

d) Optional presentation of visual, auditory, and generators, and the optional display of the multivariate 
somato-sensorial stimulation to the experimental subject distances. 

during EEG and MEG recording, carried out under the 45 * n accordance with a second aspect of the invention 
control of the central experimental program; there is provided a system for the three dimensional 

e) Optional recording and identification of responses tomography of activity and connectivity of brain elec- 
produced by the experimental subject during EEG and tromagnetic waves generators, including: 

MEG recording, for the inclusion of fiducial markers in a) A set of electrodes and magnetic sensors adapted to 
the recording, and for the modification of the central 50 be attached or approximated to the scalp of an experi- 
experimental program; mental subject for the detection of brain electromag- 

f) Optional real-time detection of spontaneous events netic physiological activity in the form of electroen- 
in the EEG and MEG produced by the experimental cephalogram (EEG) and magnetoencephalogram 
subject during recording, for the inclusion of fiducial (MEG) electromagnetic signals, and means for measur- 
markers in the recording, and for the modification of 55 ing the exact positions of the electrodes and sensors 
the experimental program; with respect to a reference coordinate system deter- 

g) Determination of a parametric description for the mined by certain anatomical landmarks of the subject's 
anatomy of the experimental subject's head (parametric head; 

geometry), by means of: i) exact computations based on b) Means for the amplification electromagnetic sig- 

anatomical or functional image processing of the sub* 60 nals detected at each electrode and sensor; 

ject's head, or ii) approximate computations based on a c) Means for obtaining on-line digital spatio-temporal 

small set of anatomical measurements and comparison signals consisting of EEG and MEG signals; 

with a data base of normal and abnormal variability; d) Means for the presentation of visual, auditory, and 

h) Using the parametric geometry for constructing a somato-sensorial stimulation to the experimental subject 
head phantom with all the volume conductor properties 65 during EEG and MEG recording; 

of the real head; e) Means for recording the vocal or movement re- 

i) Performing EEG and MEG measurements on the sponses produced by the experimental subject during 
head phantom due to known current dipoles located in EEG and MEG recording; 
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A central digital computer subsystem, consisting of a recording, the experimental subject can optionally be 

multitasking processor or a set of distributed processors, presented with visual, auditory, and somato-sensorial 

that includes: stimulation (4), with fiducial markers being added to the 

Means for reading the experimental subject's image recording in the usual manner for later identification, 
data in ithe form of CAT scan images, NMR images, or 5 Stimulation can be presented in the form of video- 

m the form of a small set of anatomical measurements, At the same time, responses produced by the 

and means for computmg and storing the descnptjve experimental subject (5) in the form of vocalizations or 

n^,n, !JT ^ ^ aMt ° m P 1 dcconvolution M movements can ^ rccorded Md identifledi md 

^«ncV^ * ^ u , . the spontaneous events in the EEG and MEG produced 

desSv^ n^nf ^St P T baSed f ° n 2* 10 b * the «P«imental subject as he or she makL the re- 
descriptive parametric geometry, and means for the * i , t , 

implJtationof current dipolesin the corresponding *Z °C * , * ^ T*" 

neural tissue volume of therfiantom; mental plan. The events can be paroxistic events or 

Means for pro^amming and for the control ofexper- *P ha ^ythm desynchromzation events, and the events 
iments that comprise stimulation of the experimental 15 ^ either J* defmed detennimstically, statistically, or 
subject, recording of the subject's responses, detection m * erms °" fx ^ y l°? c * . , . , 
and recording of special EEG and MEG events, and . Anatomical and functional information about the 
simultaneous recording of the digitized electromagnetic bead * mcludm S aspects such as geometry of the differ- 
signals; constituent tissues (e.g., brain, skull, scalp), location 

Means for pre-processing the recorded electromag- 20 and orienta ^ on of the cortex and of other neuronal 
netic signals for artifact and noise elimination; aggregates, is summarized in what is here termed de- 

Means for estimating event related components scriptive parametric geometry (8). Each surface bound- 
(ERCs); ary can in general be expressed as an implicit function 

Means for computing the ERCs sufficient statistics; F(x,y ( z)=0, where the variables x, y, and z represent 

Means for estimating the additive non-white spatio- 25 common cartesian coordinates. The descriptive para- 
temporal noise due to diffuse generators; metric geometry of an experimental subject consists of 

Means for performing tests of hypotheses about the the finite set of coefficients {Cmijk} used in approximat- 
goodness of fit of the estimated inverse solution; ing each boundary m (m= 1,2, . . . ,N« where is the 

Means for estimating the localizations, orientations, number of boundaries) in terms of an expansion of the 
activities, and connectivities of the ERCs generators; 30 form: 

Means for comparing the ERCs generators charac- 
teristics with a normative data base and means for com- 
puting multivariate metrics; ? ? | c m gkB&)Bfy)BtiZ) 

Means for the visual display of ERCs generators ' J 

characteristics and of the multivariate metrics, 35 wbere {Bn} h a finile set of R — ► R basis functions 

BRIEF DESCRIPTION OF THE DRAWINGS (n=l,2, . . . Ny; N/is number of basis functions). 

In order that the invention be better understood, F ° r each boundarv » ^ expansion is fitted (in a least 
further detailed description follows with reference to SS« ^ to J^*° m anatomical images (7) e.g., 
accompanying drawings in which: 40 CAT ™ ages or NMR ima S«0 or from functional im- 

FIG. 1) is a schematic drawing of one embodiment of ages <«*• J^T »»?S» or SPECT of the experi- 

the present invention. mental sub J e ct's head: 

FIG. 2) illustrates an example of a visual display of 
the characteristics of diffuse and concentrated genera- _ . 2 

t0 5£ *~ m?^ ^P*™*™** **J«*- 45 Min 2 [FJX m YM -222 Q^ffC^WMd J 1 

FIG. 3) illustrates an example of visual display of the cUV i j k "** J J 

characteristics of diffuse and concentrated generators in 

a patient with a cerebrovascular accident. where u= 1, 2, . . . , M m ; and M m is the number of sam- 

DETAILED DESCRIPTION OF THF plcd P ? ktS ° n m " th boUDdarv * Alternatively, the 

P^reRRFD F^r^MF^r^ 50 ^P***™ coefficients can be statistically predicted 

PREFERRED EMBODIMENTS from a minimum set of easily measured anatomical pa- 

A plurality of sensors (1) (electrodes for EEG, mag- rameters (7) by means of the information contained in a 

netometers or gradiometers for MEG) are placed on or normative data base with a representative sample of 

in proximity to the experimental subject's scalp, for the experimental subjects covering the typical population: 
detection of brain electromagnetic signals, which are 55 

generated due to neuronal sources. The sensors are C m ,jk=Gmyk (0) 
placed in order to maximize the amount of information 

about the generators. Due to variability of head shape where Qmijk is the statistically estimated R*-* R func- 

and size in the human population, measurements of the tion, and 0c ft 5 is the anatomical measurement parame- 

exact positions of the sensors are required, with respect 60 ter vector. 

to a reference coordinate system determined by certain Important basis functions B to be used are: 

anatomical landmarks of the individual subject's head. a) The Fourier basis B„(v)=exp (ivo n ), where 

The electromagnetic signals are amplified (2) to the i=(— 1)1, and v, a>„€ R. 

dynamic range of an analog to digital converter (3), b) The Karhunen-Loeve basis for each type of bound- 

which converts the signals into numbers that are stored 65 ary surface, 

in the memory of a digital computer. c) Spatial spline bases. 

The recording of electromagnetic signals is carried The descriptive parametric geometry previously ex- 
out according to a central experimental plan. During plained can be used for computing the anatomical de- 
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convolution operator { 0, E } (B0) f which is defined in spatio-temporal properties of the stochastic ERCs in 

the following equations: terms of sufficient statistics (113). The sufficient statistics 

can consist of cumulants of any order, or multiple time 
series parametric models, either in the time, frequency, 

. _ , . , 5 or time-frequency domain. Alternatively, Karhunen- 

where Vc €P and Be €9 are the electric potential and hxvc % representations can be used for fitting the 

magnetic field component measurement vectors (the statistics m stationary or n0n . stat ionary, linear 

vector elements correspond to measurements made at or non . Imear modeK 

different sensor positions), respective y ; V ro P.P and The generator modeI for E RCs is: 
Boo"* correspond to electnc potential and magnetic 10 
field component values at the same sensor positions, and 

due to the same neuronal generators, in an infinite ho- Vt = ? + E „ M Bt _ , ^ MjfcGfa + Ejft 
mogeneous medium; and 0«R Fx ^and EcR? x fare * * 
transfer coefficient matrices, which define the anatomi- 
cal deconvolution operator. J5 where the vectors V* and B,€ 8* are the time (t) 
The anatomical deconvolution operator { D, E } can dependent electric potential and magnetic field compo- 
be computed by the method and system described in nent measurements, respectively; $*«RP x3 and 
CUBAN patent application 4/91, or alternatively by the »|f 3 denote the electric and magnetic transfer 
new method described in the present invention. The coefficient matrices for the k-th dipole source, respec- 
descriptive parametric geometry is used for construct- 2 q tively; M*c R 3 is the k-th dipole's moment; gfcr is the 
ing a head phantom (9) with all the volume conductor time varying k-th dipole's intensity; E vt € RPand Est R ? 
properties of the experimental subject's head, said vol- are noise vectors; k= 1,2, . . . ,N</and Njis the number 
ume conductor properties consisting of the geometry of dipole sources. This model has the following charac- 
and conductivity profile of the different constituent teristics: 

tissues (e.g., brain, skull, scalp). Electric potential (V) 25 a ) xh e ERCs are due to N rf dipoles, with fixed posi- 
and magnetic field component (B) measurements are t j ons ^ orientations, which may be known or un- 
performed on the head phantom due to a plurality of known, and with unknown time varying dipole intensi- 
implanted current dipoles or conductivity profiles (one ties di le mtens ities are modeled as stochastic 
at a tune) with known locations and moments (dipoles processes 

located in the corresponding neural ^tissue volume). b) Measur ements consist of the ERCs (due to dipoles 

Also, theoretical EEG (Voe) and MEG (Boo) values prcvious i y described) contaminated by additive- 

are computed inthe corresponding mfimte homogene- £ modeled as a stochastic 

cms medium^ The set of vectors {(V,B,Vcc,B«,)} not necessarily „ spatio-temporal white noise. 

t^Sofe , Particular model for tlie noise Pillows for 

cal deconvolution operator { d , E } in a least squares 35 ^ S eneral ^-colored pror^rtes^ with homogeneous 

sense . r ^ isotropic spatial properties due to diffuse generators. 

Inverse solutions (14) are computed for the source 
model in infinite homogeneous medium, based on the 

[n sufficient statistics of the measured data transformed to 

$ IK^«)f — A»0/I1 J jAnd 40 infinite homogeneous medium by means of the anatomi- 
cal deconvolution operator (12), and taking into ac- 
count the constraints imposed on the generators (11). 

Minf x || (B K )t - - EtV)i\\ 2 | As a simple particular example, consider the zero 

D L ' J mean linear stationary Gaussian case for spontaneous 

45 EEG activity, where the cross spectral matrices (fre- 
where i=l,2, . . . ,Nrfand || X || denotes the norm of the quency domain second order moments) constitute the 
vector X. sufficient statistics. Assuming that the ERCs concen- 
A further important use of the descriptive parametric trated generators and the noise diffuse generators are 
geometry is the determination of anatomical and func- statistically independent, the model cross spectral rela- 
tional constraints for the localizations, orientations, 50 trix at frequency to for the electric potential (2 V ) is: 
activities, and connectivities of the brain electromag- 
netic waves generators (HI) (generator constraints). £ v -4&& j fc£°6+x ll 
Generator constraints are necessary for obtaining a 

unique inverse solution (14). For example, if the mea- where 2.y&C p*P; fanPxQNfy <^=(£io<f>2"<fw</); 

sured EEG and MEG activity is known beforehand to 55 m £ r (3Wrf)xAtf 
be generated only by cortical sources, then the genera- 
tors can be located only on the cortical surface with 
orthogonal orientations. 

Digital pre-processing (12) of the recorded EEG and 
MEG is necessary before proceeding to obtain inverse 60 

solutions. The two main steps consist of: y 

a) Artifact and noise elimination. 

b) Separation of EEG and MEG samples related to 
fiducial markers (e.g. stimuli, subject's response, or 
spontaneous EEG and MEG events) for obtaining 65 
event related components (ERCs). 

Statistical hypothesis tests are performed for deter- 2g€ CNdx Nd\ % tne generators cross spectral matrix; and 

mining the most adequate numerical description of the 2n« R^x'is the pure real valued symmetric homogene- 
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ous isotropic noise cross spectral matrix, where, in the 
case of spherical geometry, 2n=2AjtP*, where the 
A*£0 are unknown coefficientsm, [P*]#=P*(cos (y/y)), 
P* being the Legendre polynomials of order k, and yy 
is the angle between the i-th and j-th electrodes. 5 

Inverse solutions canJbe obtained by at least squares 
criterion in which Tr(2 v — V— 2 V )* is minimized with 
respect to the dipole location parameters in matrix <fr, 
the dipole orientation parameters in matrix M, the gen- 
erators cross spectral matrix 2^, and the A cpefficients a* 10 
of the noise cross spectral matrix, where 2 v is the sample 
electric potential cross spectral matrix. Alternatively, 
inverse solutions can be obtained by maximizing the 
likelihoodfunction, which is equivalent to rninimizing 
tr<2,-i£ r )-det(2 v ->i v ). Independently of the estima- 15 
tion method used, the inverse solutions must be ob- 
tained under the generator constraints. Once the loca- 
tions and orientations of the ERCs dipoles are esti- 
mated, their activities G r *=(gi+g2+"gMft) can be esti- 
mated at each time instant as: 2Q 

which is obtained by either maximum likelihood or 
weighted least squares methods. 

Summarizing the inverse solution example previously 25 
described, the following generator characteristics were 
computed for a given number of dipoles and a given 
number of coefficients a* in the expansion of the noise 
cross spectral matrix: 

a) The activity kt for each dipole k, at each time in- 30 
stant t. 

b) The connectivities X v between all pairs of dipoles. 

c) The locations (<J>) and orientations (M) of the di- 
poles. 

d) The expansion coefficients of the homogeneous 35 
isotropic noise cross spectral matrix. 

Tests for inverse solution uniqueness (15) must be 
made before performing estimations, since under non- 
uniqueness, computations may render meaningless re- 
sults. Uniqueness may be achieved by setting more se- 40 
vere generator constraints, or by decreasing the number 
of dipoles until the model is' identifiable. 

Once the generator characteristics are estimated, a 
test for the goodness of fit of the inverse solution is 
performed (16). Also, hypotheses concerning if a given 45 
generator has significant activity, or if the connectivity 
between a pair of generators is significant, can also be 
tested (16). In all cases statistical hypothesis testing is 
based on resampling techniques such as Montecarlo, the 
Jackknife, and the Bootstrap. 50 

The estimated localizations, orientations, activities, 
and connectivities of the ERCs generators of an experi- 
mental subject are compared with those of a normal 
population (17), by means of multivariate metrics for 
measuring .distances between estimators and normative 55 
data of a sample from the normal population, taking 
properly into account the effect of covariables such as 
age in order to decrease metric variability. 

Finally, visual displays (18) are presented in the form 
of three dimensional and two dimensional images of the 60 
head, where the localizations, orientations, activities, 
and connectivities of the ERCs generators are displayed 
by coding their numerical values in terms of color, 
intensity, and graphical icons, and where optionally, the 
multivariate metrics corresponding to comparison with 65 
norms can also be displayed by superposition. 

FIG. 2 illustrates a visual display of an inverse solu- 
tion obtained from a normal experimental subject, based 
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on a spontaneous EEG recording under eyes closed, 
awake, conditions. The sufficient statistics used was the 
cross spectral matrix at the alpha peak maximum 
(w=9.75 Hz). Two dipoles were fitted, together with 
additive uncorrelated homogeneous isotropic spatial 
noise. Generators are represented as arrows in the three 
head views (back, top, and left views). Generator local- 
ization and moments are given in cartesian coordinates 
referred to a unit radius sphere, with z axis coming out 
through the vertex, x axis coming out through nasion, 
and y axis coming out through the left ear (T3 electrode 
position), connectivities are given in terms of the gener- 
ator correlation matrix. Noise characteristics are illus- 
trated as "BASE ACTIVITY", giving the values of the 
expansion coefficients of the homogeneous isotropic 
process. Note the alpha rhythm generators located in 
the occipital cortex. 

FIG. 3 illustrates a visual display of an inverse solu- 
tion obtained from an experimental subject with a later- 
alized right cerebrovascular accident, based on a spon- 
taneous EEG recording under eyes closed, awake, con- 
ditions. Analysis procedures were the same to those 
used for the normal subject of FIG. 2. Note that only 
one alpha generator lies in the normal position (left 
occipital cortex). 

We claim: 

1. A method for three dimensional tomography of 
activity and connectivity of brain electromagnetic 
waves generators, said method including the steps of: 

a) positioning a set of electrodes and magnetic sensors 
to detect and record brain electromagnetic signals 
representing physiological activity in the form of 
an electroencephalogram (EEG) and a magnetoen- 
cephalogram (MEG), and measuring exact posi- 
tions of the electrodes and sensors with respect to a 
reference coordinate system determined by certain 
anatomical landmarks of a head of an experimental 
subject; 

b) amplifying the electromagnetic signals detected at 
each electrode and sensor by connecting an ampli- 
fier to each electrode and sensor; 

c) obtaining on-line digital spatio-temporal signals, 
consisting of said EEG and MEG, by connecting 
analog-digital converters to each amplifier, and 
digitizing all data as it is gathered; 

d) determining a parametric description for an anat- 
omy of the experimental subject's head to obtain a 
descriptive parametric geometry; 

e) using said descriptive parametric geometry for 
constructing a head phantom with volume conduc- 
tor properties of the experimental subject's head; 

f) performing EEG and MEG measurements on said 
head phantom due to known current dipoles lo- 
cated in a corresponding neural tissue volume, for 
determining a linear operator which transforms 
original EEG and MEG measurements into equiv- 
alent infinite homogeneous medium measurements 
(anatomical deconvolution); 

g) using said descriptive parametric geometry for 
determining anatomical and functional constraints 
for localizations, orientations, activities, and con- 
nectivities of the brain electromagnetic waves gen- 
erators (generator constraints); 

h) digitally pre-processing the EEG and MEG for 
artifact and noise elimination, and for separation of 
EEG and MEG samples related to said fiducial 
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markers, for obtaining event related components 
(ERCs); 

i) statistically analyzing said ERCs for determining 
the most adequate numerical description of spatio- 
temporal properties in terms of sufficient statistics; S 

j) computing activities and connectivities of the 
ERCs generators, based on a static solution to an 
inverse electromagnetic problem, under said gener- 
ator constraints, using sufficient statistics for the 
ERCs transformed to an infinite homogeneous 10 
medium by means of said anatomical deconvolu- 
tion; 

k) in case said generator constraints do not allow a 
unique solution to the inverse problem, decreasing 
a number of ERCs generators sufficiently to allow IS 
for proper identifiability of the inverse problem; 

1) statistically evaluating a goodness of fit of the static 
solution, taking into account an existence of col- 
ored spatial and temporal noise, and including sta- 
tistical hypotheses for testing an absence of activity 20 
and connectivity of the ERCs generators; and 

m) visually displaying three dimensional and two 
dimensional images corresponding to the localiza- 
tions, orientations, activities, and connectivities of 
the ERCs generators. 25 

2. A method as claimed in claim 1 wherein events are 
detected in the EEG and MEG for recording and modi- 
fication of an experimental sequence. 

3. A method as claimed in claim 2 wherein said events 
are defined detenninistically. 30 

4. A method as claimed in claim 2 wherein said events 
are defined statistically. 

5. A method as claimed in claim 2 wherein said events 
are determined in terms of fuzzy logic. 

6. A method as claimed in claim 2 wherein said events 35 
are paroxistic events. 

7. A method as claimed in claim 2 wherein said events 
are alpharhythmic desynchronization events. 

8. A method as claimed in claim 1 wherein series 
expansions are established in terms of three dimensional 40 
basis functions for the description of contours, physical 
properties, and metabolic properties of the head, pro- 
ducing a descriptive parametric geometry consisting of 
sets of expansion coefficients for the basis function. 

9. A method as claimed in claim 1 wherein the de- 45 
scriptive parametric geometry establishes a priori con- 
straints on said localizations and orientations of the 
ERCs generators. 

10. A method as claimed in claim 1 wherein the de- 
scriptive parametric geometry is used for constructing 50 
the head phantom with all of said volume conductor 
properties of the experimental subject's head. 

11. A method as claimed in claim 1 wherein EEG and 
MEG measurements are performed on the head phan- 
tom, using established current dipoles located in the 55 
corresponding neural tissue volume, with known loca- 
tions, orientations and activities, where theoretical 
EEG and MEG values are computed in the correspond- 
ing infinite homogeneous medium, and where the linear 
operator in matrix form which transforms original EEG 60 
and MEG measurements into equivalent infinite homo- 
geneous medium measurements is computed in a least 
squares sense, thus defining an anatomical deconvolu- 
tion operation. 

12. A method as claimed in claim 1 wherein event 65 
related samples (ERSs) of said EEG and MEG are 
separated according to different types of fiducial mark- 
ers, and wherein said ERSs contain ERCs, which are 
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estimated under the assumption that a set of generators 
of the ERCs of said ERSs have fixed localizations and 
orientations, and stochastic activities and connectivities. 

13. A method as claimed in claim 1 wherein statistical 
hypothesis tests are used for determining a set of suffi- 
cient statistics for the description of the ERCs, where 
said set of sufficient statistics consist of cumulants of 
any order, in any domain. 

14. A method as claimed in claim 1 wherein the ERSs 
are modeled as containing additive noise due to the 
activity of diffuse generators independent of ERCs 
concentrated generators, said diffuse generators having 
the properties of a stochastic process. 

15. A method as claimed in claim 1 wherein the in- 
verse problem is solved for the activities and connec- 
tivities of neuronal generators, based on all the informa- 
tion provided by ERCs sufficient statistics, under the 
generator constraints on localizations and orientations 
provided by the descriptive parametric geometry, for 
infinite homogeneous medium measurements obtained 
by applying an anatomical deconvolution operator. 

16. A method as claimed in claim 1 wherein the good- 
ness of fit of the estimated inverse solution and tests of 
hypotheses on no activity of the ERCs generators ex- 
amined based on statistical resampling techniques. 

17. A method as claimed in claim 1 wherein the esti- 
mated localizations, orientations, activities, and connec- 
tivities of the ERCs generators of an experimental sub- 
ject are compared with those of a normal population by 
means of multivariate metrics for measuring distances 
between estimators and normative data of a sample 
from the normal population, taking into account an 
effect of convariables in order to decrease metric vari- 
ability. 

18. A method as claimed in claim 1 wherein visual 
displays are presented in the form of three dimensional 
and two dimensional images of the head, where the 
localizations, orientations, activities, and connectivities 
of the ERCs generators are displayed by coding their 
numerical values in terms of color, intensity, and graph- 
ical icons. 

19. A method as claimed in 18 wherein multivariate 
metrics corresponding to comparison with norms are 
also displayed by superposition. 

20. A method as claimed in claim 1 further compris- 
ing the step of pre visual, auditory, and somato-sensorial 
stimulation to the experimental subject during said re- 
cording of the EEG and MEG, carried out under con- 
trol of said central experimental program. 

21. A method as claimed in claim 20 further compris- 
ing the step of identifying responses produced by the 
experimental subject during said EEG and MEG re- 
cording, inclusion of fiducial markers of the type and 
time of stimulation in said recording, and modifying an 
experimental program sequence based on the responses 
produced by the experimental subject. 

22. A method as claimed in claim 21 further compris- 
ing the step of including fiducial markers in said record- 
ing for the time and type of spontaneous events detected 
in the EEG and MEG produced by the experimental 
subject during recording, and modifying the central 
experimental program based on the spontaneous events. 

23. A method as claimed in claim 1 wherein step d) 
comprises the step of determining the parametric geom- 
etry for the anatomy of the experimental subject's head 
by means of exact computations based on image pro- 
cessing of the subject's head. 
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24. A method as claimed in claim 1 wherein step d) 
comprises the step of determining the parametric de- 
scription for the anatomy of the experimental subject's 
bead based on approximate computations based on a 
small set of anatomical measurements and comparison 5 
with a data base of normal and abnormal variability. 

25. A method as claimed in claim 1 further compris- 
ing the step of computing multivariate distances be- 
tween ERG generators characteristic of said experimen- 
tal subject and of a normal population as determined by 10 
a normative data-base. 

26. A method as claimed in claim 25 further compris- 
ing the step of displaying said multivariate distances, 

27. A method as claimed in claim 1 further compris- 
ing the step of presenting stimulation sequences to the 15 
subject in the form of computerized video-games, and 
recording responses to the video-games. 

28. A method as claimed in claim 1 wherein the de- 
scriptive parametric geometry is obtained from func- 
tional images of the experimental subject's head. 20 

29. A method as claimed in claim 1 wherein the de- 
scriptive parametric geometry is obtained from a mini- 
mum set of easily measured parameters which can be 
used to statistically predict the descriptive parametric 
geometry via information contained in a normative data 25 
base with a representative sample covering a typical 
population. 

30. A method as claimed in claim 1 wherein EEG and 
MEG values are theoretically computed on a head 
phantom defined by the descriptive parametric geome- 30 
try using established current dipoles located in corre- 
sponding neural tissue volume, with known locations, 
orientations, and activities, and wherein the theoretical 
EEG and MEG values are computed in the correspond- 
ing infinite homogeneous medium, and wherein a linear 35 
operator in matrix form which transforms original EEG 
and MEG measurements into equivalent infinite homo- 
geneous medium measurements is computed in a least 
squares sense, thus defining an anatomical deconvolu- 
tion operation. 40 

31. A method as claimed in claim 1 wherein statistical 
hypothesis tests are used for determining a set of suffi- 
cient statistics for a description of stochastic event re- 
lated samples (ERSs), wherein said set of sufficient 
statistics consists of multiple time series parametric 45 
models, in a domain selected from the group consisting 
of the time, frequency, and time-frequency domains. 

32. A system for the three dimensional tomography 
of activity and connectivity of brain electromagnetic 
waves generators comprising: 50 

a) A set of electrodes and magnetic sensors adapted to 
be positioned to detect brain electromagnetic sig- 
nals representative of physiological activity in the 
form of an electroencephalogram (EEG) and a 
magnetoencephalogram (MEG), and means for 55 
measuring exact positions of the electrodes and 
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sensors with respect to a reference coordinate sys- 
tem determined by certain anatomical landmarks of 
a head of an experimental subject. 

b) Means for amplification of said electromagnetic 
signals detected at each electrode and sensor, 

c) Means for obtaining on-line digital spatio-temporal 
signals consisting of said EEG and MEG; 

d) Means for the presentation of visual, auditory, and 
somato-sensorial stimulation to the experimental 
subject during EEG and MEG recording; 

e) Means for recording vocal or movement responses 
produced by the experimental subject during EEG 
and MEG recording; 

f) A central digital computer subsystem, comprising: 
means for reading the experimental subject's image 

data, and means for computing and storing a de- 
scriptive parametric geometry, an anatomical de- 
convolution operator, and generator constraints; 

means for constructing a head phantom based on the 
descriptive parametric geometry, and means for 
the implantation of current dipoles in a correspond- 
ing neural tissue volume of the phantom; 

means controlling experiments including means for 
causing stimulation of the experimental subject, 
recording of the subject's responses, detection and 
recording of special EEG and MEG events, and 
simultaneous recording of said electromagnetic 
signals; 

means for pre-processing the recorded electromag- 
netic signals for artifact and noise elimination; 

means for estimating event related components 
(ERCs) 

means for computing sufficient statistics of the ERCs; 

means for estimating an additive non-white spatio- 
temporal noise due to diffuse generators; 

means for performing tests of hypothesis about the 
goodness of fit of an estimated inverse solution; 

means for estimating localizations, orientations, activ- 
ities, and connectivities of generators of the ERCs; 

means for comparing characteristics of the ERCs 
generators with a normative data base and means 
for computing multivariate metrics; 

means for the visual display of ERCs generators char- 
acteristics and of the multivariate metrics. 

33. A system as claimed in claim 32 wherein said 
central digital computer subsystem includes a multitask- 
ing processor. 

34. A system as claimed in claim 32 wherein said 
central digital computer subsystem includes a set of 
distributed processors. 

35. A system as claimed in claim 32 wherein said 
means for reading the experimental subject's image data 
comprises means for reading images selected from the 
group consisting of CAT scan images, NMR images, 
PET images, and anatomical measurements. 

***** 
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